Momentum-space description of three-nucleon breakup reactions including the 

Coulomb interaction 
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The Coulomb interaction between the two protons is included in the calculation of proton-deuteron 
breakup and of three-body electromagnetic disintegration of 3 He. The hadron dynamics is based on 
the purely nucleonic charge-dependent (CD) Bonn potential and its realistic extension CD Bonn + 
A to a coupled-channel two-baryon potential, allowing for single virtual A-isobar excitation. Cal- 
culations are done using integral equations in momentum space. The screening and renormalization 
approach is employed for including the Coulomb interaction. Convergence of the procedure is found 
at moderate screening radii. The reliability of the method is demonstrated. The Coulomb effect on 
breakup observables is seen at all energies in particular kinematic regimes. 

PACS numbers: 21.30.-x, 21.45.+V, 24.70.+S, 25.10,+s 



I. INTRODUCTION 

The inclusion of the Coulomb interaction in the de- 
scription of the three-nucleon continuum is one of the 
most challenging tasks in theoretical few-body nuclear 
physics . Whereas it has already been solved for elastic 
proton-deuteron (pd) scattering with realistic hadronic 
interactions using various procedures Jl|,|3,y,|j,|5|, there 
are only very few attempts 0, IS U to calculate pd 
breakup, and none of them uses a complete treatment of 
the Coulomb interaction and realistic hadronic potentials 
allowing for a stringent comparison with the experimen- 
tal data. 

Recently in Ref. [f| we included the Coulomb inter- 
action between the protons in the description of three- 
nucleon reactions with two-body initial and final states. 
The description is based on the Alt-Grassberger-Sandhas 
(AGS) equation in momentum space. The Coulomb 
potential is screened and the resulting scattering ampli- 
tudes are corrected by the renormalization technique of 
Refs. [3 El 

to recover the unscreened limit. The treat- 
ment is applicable to any two-nuclcon potential without 
separable expansion. Reference Q and this paper use the 
purely nucleonic charge-dependent (CD) Bonn potential 
|12| and its coupled-channel extension CD Bonn + A ^| , 
allowing for a single virtual A-isobar excitation and fitted 
to the experimental data with the same degree of accu- 
racy as CD Bonn itself. In the three-nucleon system the 
A isobar mediates an effective three-nucleon force and ef- 
fective two- and three-nucleon currents, both consistent 
with the underlying two-nucleon force. The treatment of 
Ref. [f| is technically highly successful, but still limited 
to the description of proton-deuteron (pd) elastic scatter- 
ing and of electromagnetic (e.m.) reactions involving 3 He 



with pd initial or final states only. This paper extends the 
treatment of Coulomb to breakup in pd scattering and to 
e.m. three-body disintegration of 3 He. In that extension 
we follow the ideas of Refs. [EEIEIj but avoid approx- 
imations on the hadronic potential and in the treatment 
of screened Coulomb. Thus, our three-particle equa- 
tions, including the screened Coulomb potential, are com- 
pletely different from the quasiparticle equations solved 
in Ref. [6j where the two-nucleon screened Coulomb tran- 
sition matrix is approximated by the screened Coulomb 
potential. In Ref. jlj] we presented for the first time 
a limited set of results for pd breakup using the same 
technical developments we explain here in greater detail. 

We have to recall that the screened Coulomb potential 
wr we work with is particular. It is screened around the 
separation r = R between two charged baryons and in 
configuration space is given by 



wr{t) — w(r) e 



-{r/Rf 



(1) 
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with the true Coulomb potential w(r) = a e /r, a e being 
the fine structure constant and n controlling the smooth- 
ness of the screening. We prefer to work with a sharper 
screening than the Yukawa screening (n — 1) of Ref. [|J. 
We want to ensure that the screened Coulomb poten- 
tial wr approximates well the true Coulomb one w for 
distances r < R and simultaneously vanishes rapidly for 
r > R, providing a comparatively fast convergence of 
the partial- wave expansion. In contrast, the sharp cut- 
off (n — ► oo) yields an unpleasant oscillatory behavior 
in the momentum-space representation, leading to con- 
vergence problems. We find the values 3 < n < 6 to 
provide a sufficiently smooth, but at the same time a 
sufficiently rapid screening around r — R like in Ref. 0; 
n = 4 is our choice for the results of this paper. The 
screening radius R is chosen much larger than the range 
of the strong interaction which is of the order of the pion 
wavelength h/m^c « 1.4 fm. Nevertheless, the screened 
Coulomb potential wr is of short range in the sense of 
scattering theory. Standard scattering theory is therefore 
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applicable. A reliable technique for solving the AGS 
equation with short-range interactions is extended in 
Ref. p| to include the screened Coulomb potential be- 
tween the charged baryons. However, the partial- wave 
expansion of the pair interaction requires much higher 
angular momenta than the one of the strong two-nucleon 
potential alone. 

The screening radius R will always remain very small 
compared with nuclear screening distances of atomic 
scale, i.e., 10 5 fm. Thus, the employed screened Coulomb 
potential wr is unable to simulate properly the physics 
of nuclear screening and, even more, all features of the 
true Coulomb potential. Thus, the approximate breakup 
calculations with screened Coulomb wr have to be cor- 
rected for their shortcomings in a controlled way. Ref- 
erences |lfl 111) give the prescription for the correction 
procedure which we follow here for breakup as we did pre- 
viously for elastic scattering, and that involves the renor- 
malization of the on-shell amplitudes in order to get the 
proper unscreened Coulomb limit. After the indicated 
corrections, the predictions for breakup observables have 
to show independence from the choice of the screening 
radius R, provided it is chosen sufficiently large. That 
convergence will be the internal criterion for the reliabil- 
ity of our Coulomb treatment. 

Configuration space treatments of Coulomb 0, 01 may 
provide a viable alternative to the integral equation 
approach in momentum space on which this paper is 
based. References 0, 13 have provided first results for 
pd breakup, but they still involve approximations in the 
treatment of Coulomb and the employed hadronic dy- 
namics is not realistic. Thus, a benchmark comparison 
between our breakup results and corresponding configu- 
ration space results is, in contrast to pd elastic scattering 
(l(j|. not possible yet. With respect to the reliability of 
our Coulomb treatment for breakup, we rely solely on 
our internal criterion, i.e., the convergence of breakup 
observables with the screening radius R; however, that 
criterion was absolutely reliable for pd elastic scattering 
and related e.m. reactions. 

Section [H] develops the technical apparatus underlying 
the calculations. Section Hill presents some characteristic 
effects of Coulomb in three-nucleon breakup reactions. 
Section llVI gives our conclusions. 



II. TREATMENT OF COULOMB 
INTERACTION BETWEEN PROTONS IN 
BREAKUP 

This section carries over the treatment of the Coulomb 
interaction given in Ref. 0] for pd elastic scattering and 
corresponding e.m. reactions, to pd breakup and to e.m. 
three-body disintegration of 3 He. It establishes a theo- 
retical procedure leading to a calculational scheme. The 
discussions of hadronic and e.m. reactions are done sep- 
arately. 



A. Theoretical framework for the description of 
proton-deuteron breakup with Coulomb 

This section focuses on pd breakup. However, the tran- 
sition matrices for elastic scattering and breakup are so 
closely connected that certain relations between scatter- 
ing operators already developed in Ref. 0] have to be 
recalled to make this paper selfcontained. 

Each pair of nucleons (/?7) interacts through the strong 
coupled-channel potential v a and the Coulomb potential 
w a . We assume that w a acts formally between all pairs 
(f3"f) of particles, but it is nonzero only for states with 
two-charged baryons, i.e., pp and pA + states. We intro- 
duce the full resolvent G^ R \Z) for the auxiliary situation 
in which the Coulomb potential w a is screened with a 
screening radius R, w a being replaced by w a R, 

GW(Z) = (Z-H -Y,v<r-J2 w °*y 1 > ( 2 ) 

where Hq is the three-particle kinetic energy operator. 
The full resolvent yields the full pd scattering state when 
acting on the initial channel state \4> a {<\i)v a i) of relative 
pd momentum q^, energy E a (qi) and additional discrete 
quantum numbers v ai and taking the appropriate limit 
Z = E a (qi) + iO. The full pd scattering state has, above 
breakup threshold, components corresponding to the fi- 
nal breakup channel states |0o(p/q/)&tv)j P/ an d q/ 
being three-nucleon Jacobi momenta, E$(pfqf) its en- 
ergy, and UQ f additional discrete quantum numbers. The 
full resolvent therefore also yields the desired S matrix 
for breakup. The full resolvent G^ (Z) depends on the 
screening radius R for the Coulomb potential and that 
dependence is notationally indicated; the same will be 
done for operators related to G^ R '(Z). Following stan- 
dard AGS notation of three-particle scattering, the 
full resolvent G^ (Z) may be decomposed into channel 
resolvents and free resolvent 

G<£\Z) = (Z - H - v a - w aR y\ (3a) 
G (Z) = (Z-Ho)- 1 , (3b) 

together with the full multichannel three-particle transi- 
tion matrices U^(Z) for elastic scattering and Uq R '(Z) 
for breakup according to 

GW(Z) = 8 Pa Gf\Z) + Gf\z)U${Z)GW{Z), 

(4a) 

G^\Z)=G (Z)U^\Z)G^\Z). (4b) 

The full multichannel transition matrices satisfy the AGS 
equations t S| 

U^\z) = S 0a Go\Z) + Y,h^ R Hz)Go(Z)U^(Z), 

(5a) 

u£\Z) = G- a \Z) +J2tW{Z)G (Z)UW(Z), (5b) 
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with 5p a = 1 — 8/3 a ; the two-particle transition matrix 
Ti R) (Z) is derived fr om the full channel interaction v n + 
w aR including screened Coulomb, i.e., 

Ti R \Z) = (v a + w aR ) + (v a + w aR )G (Z)T^(Z). 

(5c) 

In pd elastic scattering, an alternative decomposition 
of the full resolvent is found conceptually more reveal- 
ing. Instead of correlating the plane-wave channel state 
|(/> Q (q)^ Q ) in a single step to the full scattering state by 
G^ R > (Z), it may be correlated first to a screened Coulomb 
state of proton and deuteron by the screened Coulomb 
potential W£S" between a proton and the center of mass 
(cm.) of the remaining neutron- proton (np) pair in chan- 
nel a through 

G aR (Z) = {Z-H -v a - w aR - WZg)- 1 , (6a) 
G aR (Z) = Gi R Hz) + Gl R Hz)T^-(Z)Gi R \z), (6b) 
T^(Z) = W" R + W^GM(Z)T^(Z), (6c) 

where, in each channel a, w aR and W^ R are never si- 
multaneously present: When a corresponds to a pp pair, 
w aR is present and W^ R = 0; when a denotes an np pair, 
w aR — and is present. The same Coulomb cor- 

relation is done explicitly in both initial and final states. 
Thus, the full resolvent can be decomposed, in alterna- 
tive to Eq. ijjjp , as 

G^(Z) = 6p a G aR (Z) + Gp R (Z)U$(Z)G aR (Z), (7) 

yielding a new form for the full multichannel transition 
matrix 

U (R) (Z) = 5 Pa T™-(Z) + [1 + T^{Z)Gf\z)) 
xu£\z)[l + Gi R \Z)T^(Z)}. 

The reduced operator ujfa{Z) may be calculated 
through the integral equation 

tiff(Z) = S 0a [G- R {Z) + v a ] + 8p a W aR 

+ ]T(^ CT + 8 Pa Wp R )G aR {Z)Ui R \Z), 

(8b) 

which is driven by the strong potential v a and the po- 
tential of three- nucleon nature W aR = ^2 a (S a a-w aR — 
daaWaii)- This potential W aR accounts for the differ- 
ence between the direct pp Coulomb interaction and the 
one that takes place between the proton and the cm. of 
the remaining bound as well as unbound np pair. When 
calculated between on-shell screened pd Coulomb states, 
U { ^{Z) is of short -range, even in the infinite R limit. 

In the same spirit, the final breakup state to be ana- 
lyzed may not be reached in a single step; instead it may 
be correlated first to a screened Coulomb state between 



the charged particles whose corresponding Coulomb re- 
solvent keeps only the screened Coulomb interaction, 

G R (Z) = {Z -Ho-Y^w^r)' 1 - (9a) 

In the system of two protons and one neutron only the 
channel a = p, corresponding to a correlated pp pair, 
contributes to G R (Z), 

G R (Z) = G (Z) + Go(Z)T pR (Z)G (Z), (9b) 
T pR (Z) = w pR + w pR Gq(Z)T pR (Z), (9c) 

making channel p the most convenient choice for the de- 
scription of the final breakup state. Thus, for the pur- 
pose of pd breakup, a decomposition of the full resolvent, 
alternative to Eq. (|4b[) is 

G( R \Z) = G R {Z)U { R \z)G aR {Z), (10a) 
GW{Z) = G (Z)[1 + T pR (Z)G (Z)]U^(Z) 

x[1 + GW(Z)T^(Z)]GW(Z), (10b) 

where the full breakup transition matrix may be written 
as 

U$(Z) = [1 + T PR (Z)G (Z)}U^\Z) 
x[l + G[ R Hz)T^(Z)}, 

The reduced operator U^ R \z) may be calculated 
through quadrature 

U£\Z) = G- R {Z) + v a + J2 v a G aR {Z)u£){Z) 

(lib) 

from the correspondingly reduced operator U^{Z) of 
elastic scattering. In the form i|lla|) for the full breakup 
transition matrix the external distortions due to screened 
Coulomb in the initial and final states are made explicit. 
On-shell the reduced operator Uq R \z) calculated be- 
tween screened Coulomb distorted initial and final states 
is of finite range and has two contributions with slightly 
different range properties: 

(a) The contribution G~ R (Z) + v a , when calculated 
on-shell between initial pd and final three-nucleon states, 
becomes the three-nucleon potential W aR and is the 
longest-range part of breakup, since the np pair is cor- 
related by the hadronic interaction only in the initial pd 
state. The corresponding contribution in Ref. Q is called 
the pure Coulomb breakup term. 

(b) The remaining part v a G aR {Z)l]^ {Z) is of 
shorter range, comparable to the one of the reduced op- 
erator Up R \z) for elastic pd scattering. 

In the full breakup operator U^ R \z) the external 
distortions show up in screened Coulomb waves gener- 
ated by [1 + G {R} (Z)T^ (Z)] in the initial state and by 
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[1 + T p ji(Z)Go(Z)] in the final state; both wave func- 
tions do not have proper limits as R — > oo. Therefore 

( R) 

Uq u (Z) has to get renormalized as the corresponding 
amplitude for pd elastic scattering 0, , in order to ob- 
tain the results appropriate for the unscreened Coulomb 
limit. According to Refs. the full breakup tran- 

sition amplitude for initial and final states \<p a ( ( h) l 'ai) 
and |0 o (p/q/)^o / ), E a (qi) = E Q (p f q s ), referring to the 
strong potential v a and the unscreened Coulomb poten- 
tial w a , is obtained via the renormalization of the on- 

( R) 

shell breakup transition matrix U^ a (E a (qi) + iO) in the 
infinite R limit 

= J im {^ 2 b/)(</'o(p/q/)^o / l (12) 

it— »oo v ' 

X U^ ) (E a (q l )+tO)\M^> at )Z R hq 1 )}, 

where Zn(qi) and za(pf) are pd and pp renormalization 
factors defined below. 

As in Ref. we choose an isospin description for the 
three baryons in which the nucleons are considered identi- 
cal. The two-baryon transition matrix (Z) becomes 
an operator coupling total isospin T = ^ and T = § 
states as described in detail in Ref. Instead of the 
breakup amplitude given by Eq. l|12fl we have to use the 
properly symmetrized form 

(<^0 (P/q/ ) VQ f I U 1 <f> a (qi)v ai ) 

= ^Mp/q/K|£W|0a(q*Hi): (13a) 

<J 

(<Mp/q/)^o/l^o|0 

_1 

= jI^^b 2 (Pf)(MPfClfH f \ (13b) 
x U^ R \E a { qi ) +iQ)\4> a {c h )v az )Z^{q l )} 

with U^ R) (Z) = U$\Z) + U$(Z)P 231 + U^ R) (Z)P 312 
for the calculation of observables, (aflj) being cyclic and 
P231 and P312 being the two cyclic permutations of (a/37). 

The symmetrized breakup transition matrix Uq R \z) fol- 
lows by quadrature 

U (R \z) = {l + P)G^(Z) 

+ (l + P)Ti R \Z)G (Z)U^(Z) 

from the symmetrized multichannel transition matrix 
UW(Z) = uL R J(Z) + uif (Z)P 231 + u£\z)P 312 of elas- 
tic pd scattering, satisfying the standard symmetrized 
form of the AGS integral equation (|5ap . i.e., 

U {R \Z) = PG^(Z) + PT^(Z)G (Z)U iR \z), 

(14b) 

With P = P 2 31 + p312- 

The renormalization factors Z^qi) and zn(pf) in the 
initial and final channels are diverging phase factors de- 
fined in Ref. [lfj for a general screening and calculated 



in Refs. 0, 0] for the screened Coulomb potential of 
Eq. QJ, i.e., 

z R { qi ) = e**)-cy»] t (15a) 

«fl(P/) - e- 2 ^^ & t R )~ c M 7 (15b) 

K(qi) = a e M/qi and K(pf) = a e /i/pf being the pd and 
pp Coulomb parameters, M and [i the reduced pd and 
pp mass, C ss 0.5772156649 the Euler number, and n the 
exponent in Eq. Q). In pd elastic scattering, the renor- 
malization factors were used in a partial-wave dependent 
form which yielded a slight advantage on convergence 
with R compared to the partial-wave independent form 
ijTSJl. In breakup, the operator ri R] \Z)G a {Z)U {W > (Z) in 
Eq. i|14a|) is calculated in a partial- wave basis, but the on- 

( R) 

shell elements of the full breakup operator Uq (Z) are 
calculated in a plane wave basis. Therefore the renormal- 
ization is only applicable in the partial- wave independent 
form of Eq. ljTo)l, 

The limit in Eq. (|13b|) has to be performed numeri- 
cally, but, due to the finite-range nature of the breakup 
operator, the infinite R limit is reached with sufficient 
accuracy at rather modest screening radii R. Further- 
more, the longer-range pure Coulomb breakup part which 
after symmetrization reads [1 -I- T p n(Z)Go(Z)]Pv a [l + 

Ga (Z)T^(Z)] and the remaining shorter-range part 
can be renormalized with different screening radii, since 
the limit in Eq. I|13b|l exist for them separately. The 
limit for the pure Coulomb breakup part can even 
be carried out explicitly, since the renormalization of 
the screened Coulomb waves yields the corresponding 
unscreened Coulomb waves accessible in configuration 
space; thus, the integral can be carried out numerically 
in configuration space as was done indeed in Ref. 
However, we find such a procedure unnecessary when 
our standard screening function is used. In fact, in 
most cases there is even no necessity for splitting the 
full breakup amplitude into pure Coulomb and Coulomb 
modified short-range parts, the only exception being the 
kinematical situations characterized by small momentum 
transfer in the pp subsystem which are sensitive to the 
Coulomb interaction at larger distances. 

The practical implementation of the outlined calcu- 
lational scheme faces a technical difficulty. We solve 
Eq. jHEIl in a partial- wave basis. The partial- wave ex- 
pansion of the screened Coulomb potential converges 
rather slowly. In this context, the perturbation the- 
ory for higher two-baryon partial waves developed in 
Ref. i s a very efficient and reliable technical tool for 
treating the screened Coulomb interaction in high par- 
tial waves. We vary the dividing line between partial 
waves included exactly and perturbatively in order to 
test the convergence and thereby establish the validity 
of the procedure. Furthermore, the partial-wave conver- 
gence becomes slightly faster when replacing lowest-order 
screened Coulomb contributions in Uq (Z) by the re- 
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spective plane-wave results, i.e., 1-2 



U ( H \Z) = [U^>(Z) - (1 + P)w afl P] + (1 + P)w afl P, 

(16) 

where the first term converges with respect to partial 
waves faster than L/q (Z) itself and the second term is 
calculated without partial-wave decomposition. 

With respect to the partial-wave expansion in the ac- 
tual calculations of this paper, we obtain fully converged 
results by taking into account the screened Coulomb in- 
teraction in two-baryon partial waves with pair orbital 
angular momentum L < 15; orbital angular momenta 
9 < L < 15 can safely be treated perturbatively. The 
above values refer to the screening radius R = 30 fm; for 
smaller screening radii the convergence in orbital angular 
momentum is faster. The hadronic interaction is taken 
into account in two-baryon partial waves with total an- 
gular momentum I < 5. Both three-baryon total isospin 
T = I an( i T = f states are included. The maximal 
three-baryon total angular momentum J considered is 

61 
2 ■ 

FiguresQ-Elstudy the convergence of our method with 
increasing screening radius R according to Eq. (|13b|l . All 
the calculations of this section are based on CD Bonn as 
the hadronic interaction. The kincmatical final-state con- 
figurations are characterized in a standard way by the po- 
lar angles of the two protons and by the azimuthal angle 
between them, (01,02, ^2 — (fi). We show several char- 
acteristic examples referring to pd breakup at 13 MeV 
proton lab energy and at 130 MeV deuteron lab energy. 
The convergence is impressive for the spin-averaged dif- 
ferential cross section as well as for the spin observables 
in most kinematical situations as demonstrated in Figs.^ 
and [21 The screening radius R = 20 fm is sufficient; 
only in the top plot of Fig. the curves for R = 20 fm 
and R = 30 fm are graphically distinguishable. The ex- 
ception requiring larger screening radii is the differential 
cross section in kinematical situations characterized by 
very low pp relative energy E pp , i.e., close to the pp final- 
state interaction (pp-FSI) regime, as shown in Fig. [3J In 
there, the pp repulsion is responsible for decreasing the 
cross section, converting the pp-FSI peak obtained in the 
absence of Coulomb into a minimum with zero cross sec- 
tion at pf = 0, i.e., for E pp = 0. A similar convergence 
problem also takes place in pp scattering at very low en- 
ergies as discussed in Ref. 5]. In fact, screening and 
renormalization procedure cannot be applied at pf = 0, 
since the renormalization factor Zr(pj — 0) is ill-defined. 
Therefore an extrapolation has to be used to calculate 
the observables at pf — 0, which works pretty well since 
the observables vary smoothly with. pf. In Fig.[3]thc fully 
converged result would start at zero for E pp = 0. 

The seen Coulomb effects and their physics implica- 
tions are discussed in Sec. IIIII 
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FIG. 1: Convergence of the pd breakup observables with 
screening radius R. The differential cross section and the 
proton analyzing power A y (N) for pd breakup at 13 MeV 
proton lab energy are shown as functions of the arclength S 
along the kinematical curve. Results for CD Bonn potential 
obtained with screening radius R = 10 fm (dotted curves), 
20 fm (dash-dotted curves), and 30 fm (solid curves) are com- 
pared. Results without Coulomb (dashed curves) are given as 
reference for the size of the Coulomb effect. 



B. Three-body e.m. disintegration of 3 He 

For the description of the considered e.m. processes 
the matrix element (ip^ (P/Q/)^/ b 7 '(Q, K + )|P) of the 
e.m. current operator between the three-nucleon bound 
state and the breakup scattering state has to be calcu- 
lated. The calculation of that matrix element without 
Coulomb and the meaning of the momenta Q and K+ 
are discussed in great length in Refs. 0,H(j- This sub- 
section only discusses the modification which arises due 
to the inclusion of the Coulomb interaction between the 
charged baryons. Coulomb is included as a screened po- 
tential and the dependence of the bound and scattering 
states, i.e., \BW) and |4 ±Ki V/q/K), on the screen- 
ing radius R is notationally made explicit. In analogy to 
pd breakup, the current matrix element referring to the 
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FIG. 2: Convergence of the pd breakup observables with 
screening radius R. The differential cross section and the 
deuteron analyzing power A xx for pd breakup at 130 MeV 
are shown. Curves as in Fig. 
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FIG. 3: Convergence of the pd breakup observables with 
screening radius R. The differential cross section for pd 
breakup at 13 MeV proton lab energy in the pp-FSI config- 
uration is shown as function of the relative pp energy E vv . 
Results obtained with screening radius R = 10 fm (dotted 
curves), 20 fm (dashed-double-dotted curves), 30 fm (dashed- 
dotted curves), 40 fm (double-dashed-dotted curves), and 
60 fm (solid curves) are compared. Results without Coulomb 
(dashed curves) are given as reference for the size of the 
Coulomb effect. 

unscreened Coulomb potential is obtained via renormal- 
ization of the matrix element referring to the screened 
Coulomb potential in the infinite R limit 

<^- ) ( P/ q / KI^(Q,K + )|i?) 

= lim {z R H Pf )(4- ){R \p f c if )v af \r(Q,K + )\BW)}. 

rt — >oc 

(17) 

The renormalization factor Zu{pr) is the same as used in 
pd breakup for the final state. Due to the short-range na- 
ture of ji*(Q,K + )\BW) the limit R — > oo is reached with 
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FIG. 4: Convergence of the 3 He(7,pp)n reaction observables 
with screening radius R. The differential cross section and 
the target analyzing power A y at 15 MeV photon lab energy 
are shown. Curves as in Fig. 



sufficient accuracy at finite screening radii R. The pres- 
ence of the bound-state wave function in the matrix ele- 
ment strongly suppresses the contribution of the screened 
Coulomb interaction in high partial waves, i.e., two- 
baryon partial waves with orbital angular momentum 
L < 6 are sufficient for convergence. The other quantum- 
number related cutoffs in the partial-wave dependence of 
the matrix element are the same as in Refs. [l!l|2(|, i.e., 
I < 4, J < y for photoreactions, and / < 3, J < ^ for 
inelastic electron scattering from 3 He. All calculations 
include both total isospin T = \ and T — | states. 

Figures 0] and study the convergence of our method 
with increasing screening radius R for the three-body 
photodisintegration of 3 He at 15 MeV and 55 MeV pho- 
ton lab energy. The calculations are again based on CD 
Bonn as the hadronic interaction and the currents from 
Refs. 0,[2(|. We show the differential cross section and 
the target analyzing power A y for selected kinematic con- 
figurations. The convergence is again extremely good 
and quite comparable to pd breakup; the screening ra- 
dius R — 20 fm is fully sufficient in most cases. The 
only exceptional cases, as in the pd breakup, are pp-FSl 
regimes as shown in Fig. [S] The convergence with in- 
creasing screening radius R is the same for three-body 
electrodisintegration of 3 He; we therefore omit a corre- 
sponding figure. 
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FIG. 5: Convergence of the 3 He(7, pn)p reaction observ- 
ables with screening radius R. The differential cross section 
at 55 MeV photon lab energy in the pp-FSI configuration is 
shown. Curves as in Fig. [3] 



III. RESULTS 

We base our calculations on the two-baryon coupled- 
channel potential CD Bonn + A with and without 
Coulomb and use the CD Bonn potential with Coulomb 
as purely nucleonic reference. We use the charge and cur- 
rent operators of Refs. |19U20|. appropriate for the under- 
lying dynamics. In contrast to Ref. |5j we do not include 
one-nucleon relativistic charge corrections for photoreac- 
tions, since their effect on the considered observables is 
very small. 

Obviously, we have much more predictions than it is 
possible to show. Therefore we make a selection of the 
most interesting predictions which illustrate the message 
we believe the results tell us. The readers are welcome 
to obtain the results for their favorite data from us. 



A. Proton-deuteron breakup 

FiguresEl-|5|give our results for the fivefold differential 
cross section at 10.5 MeV, 13 MeV, 19 MeV, and 65 MeV 
proton lab energies in the standard space star, collinear, 
quasifree scattering (QFS), and np final-state interac- 
tion (np-FSI) configurations, for which there is available 
experimental data. Though the inclusion of Coulomb 
slightly improves the agreement with data in the space 
star configurations in Fig. El the Coulomb effect is far 
too small to reproduce the difference between pd and nd 
data and to resolve the so-called space star anomaly at 13 
MeV. The inclusion of Coulomb clearly improves the de- 
scription of the data around the collinear points at lower 
energies, i.e., at the minima in Fig. \7\ The remaining 
discrepancies around the peaks are probably due to the 
finite geometry, not taken into account in our calculations 
owing to the lack of information on experimental details, 
but may also be due to the underlying hadronic interac- 



tion. The inclusion of Coulomb decreases the differential 
cross section around the QFS peaks, i.e., around the cen- 
tral peaks in Fig. |SJ those changes are supported by the 
data at lower energies. In the np-FSI configurations of 
Fig-IHthe Coulomb effect is rather insignificant. 

The Coulomb effect on proton analyzing powers in the 
considered kinematical configurations is usually small on 
the scale of the experimental error bars. We therefore 
show in Fig. ^] only few collinear configurations. 

Recently, pd breakup has been measured at 130 MeV 
dcutcron lab energy in a variety of kinematical configura- 
tions j2^. In some of them we find significant Coulomb 
effects for the differential cross section as well as for 
the deuteron analyzing powers. Examples are shown 
in Figs. El an d [HI By and large the agreement be- 
tween theoretical predictions and experimental data is 
improved. The pp-FSI repulsion is responsible for lower- 
ing the peak of the differential cross section in the config- 
uration (15°, 15°, 40°) in Fig. El left, where the relative 
pp energy is rather low at the peak. In contrast, the rela- 
tive pp energy gets considerably increased as one changes 
the azimuthal angle to 160° in Fig.EJright, leading to an 
increase of the differential cross section due to Coulomb. 
Since the total breakup cross section at this energy, cor- 
responding to 65 MeV proton lab energy, is almost unaf- 
fected by Coulomb, as shown in Fig. 1131 one may expect 
in given configurations an increase of the cross section 
due to Coulomb to compensate for the sharp decrease of 
the cross section in the vicinity of pp-FSI points. Figure 
II 21 shows deuteron tensor analyzing powers A xx and Ayy 
with moderate Coulomb effect for the same configura- 
tions for which experimental data may become available 
soon |2S|. 

Compared to the results of Ref. @ based on a simple 
hadronic S'-wave potential, we see a rough qualitative 
agreement in most cases. Quantitatively, the Coulomb 
effect we observe is smaller than the one of Ref. @. 

Figures lol - 1 1 31 recall also the A-isobar effect on observ- 
ables, which, in most cases we studied, is much smaller 
than the Coulomb effect. As expected, the A-isobar ef- 
fect on polarization observables is more significant than 
on the differential cross sections, which confirms previous 
findings [l3[ . 



B. Three-body e.m. disintegration of 3 He 

Experimental data for three-body photodisintegration 
of 3 He are very scarce; we therefore show in Fig. 1141 only 
two examples referring to the semiinclusive 3 He(7, pn)p 
reaction at 55 MeV and 85 MeV photon lab energy. 
The semiinclusive fourfold differential cross section is ob- 
tained from the standard fivefold differential cross section 
by integrating over the kinematical curve S. For scat- 
tering angles corresponding to the peak of the fourfold 
differential cross section the region of the phase space 
to be integrated over contains pp-FSI regime where the 
pp-FSI peak obtained without Coulomb is converted into 
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FIG. 6: Differential cross section for space star configurations 
as function of the arclength S along the kinematical curve. 
Results including A-isobar excitation and the Coulomb inter- 
action (solid curves) are compared to results without Coulomb 
(dashed curves). In order to appreciate the size of the A- 
isobar effect, the purely nucleonic results including Coulomb 
are also shown (dotted curves). The experimental pd data 
(circles) are from Ref. [gj at 10.5 MeV, from Ref. |23 at 13 
MeV, from Ref. |23| at 19 MeV, from Ref. Igl at 65 MeV, 
and nd data at 13 MeV (squares) from Ref. |25| . 



FIG. 7: Differential cross section for collinear configurations. 
Curves and experimental data as in Fig.|S] except for 65 MeV 
data from Ref. |2jj. 



a minimum as shown in Fig. [3] Therefore the fourfold 
differential cross section in Fig. [21 is also significantly 
reduced by the inclusion of Coulomb, clearly improving 
the agreement with the data. A similar Coulomb effect of 
the same origin is shown in Fig. 1151 for the semiinclusive 
threefold differential cross section for 3 He(7, n)pp reac- 
tion at 15 MeV photon lab energy. In contrast, the pho- 
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FIG. 8: Differential cross section for QFS configurations. 
Curves and experimental data as in Fig.|H] except for 65 MeV 
data from Ref. 
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FIG. 9: Differential cross section for np-FSI configurations. 
Curves and experimental data as in Fig. |H] 



we show 3 He inclusive longitudinal and transverse re- 
sponse functions Rl and Rt as examples. Though the 
Coulomb effect may be large in particular kinematic re- 
gions, it is rather insignificant for the total cross section 
and therefore also for response functions. Only the trans- 
verse response function near threshold is affected more 
visibly as shown in Fig. I17I at higher momentum trans- 
fer there is also quite a large A-isobar effect. 



IV. SUMMARY 



ton analyzing power remains almost unchanged by the 
inclusion of Coulomb. The experiment measuring this 
reaction is in progress [3lJ > but the data are not available 
yet. Again the importance of A-isobar degree of freedom 
is considerably smaller than the effect of Coulomb. 

The available data of three-nucleon electrodisintegra- 
tion of 3 He refer to fully inclusive observables. In Fig.llbl 



In this paper we show how the Coulomb interac- 
tion between the charged baryons can be included 
into the momentum-space description of proton-deuteron 
breakup and of three-body e.m. disintegration of 3 He 
using the screening and renormalization approach. The 
theoretical framework is the AGS integral equation 
The calculations are done on the same level of accu- 
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FIG. 10: Proton analyzing power for collinear configurations 
at 13 MeV and 65 MeV proton lab energy. Curves and ex- 
perimental data as in Fig. 
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FIG. 11: Differential cross section for pd breakup at 130 MeV 
deuteron lab energy. Curves as in Fig. |S| The experimental 
data are from Ref. pflj . 
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FIG. 12: Deuteron analyzing powers for pd breakup at 130 
MeV deuteron lab energy. Curves as in Fig. |S| 
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FIG. 13: Total cross section for pd breakup as function of 
the proton lab energy Curves as in Fig. [5] The experimental 
data are from Ref. 13011 . 



racy and sophistication as for the corresponding neutron- 
deuteron and 3 H reactions. The conclusions of the paper 
refer to the developed technique and to the physics re- 
sults obtained with that technique. 

Technically, the idea of screening and renormalization 
is the one of Refs. |6l lloL[rT| . However, our practical real- 
ization differs quite significantly from the one of Ref. |(| : 

(1) We use modern hadronic interactions, CD Bonn 
and CD Bonn + A, in contrast to the simple iS-wave 
separable potentials of Ref. Our use of the full po- 
tential requires the standard form of the three-particle 
equations, different from the quasiparticle approach of 
Ref. @. 

(2) We do not approximate the screened Coulomb tran- 
sition matrix by the screened Coulomb potential. 

(3) The quasiparticle approach of Ref. |(j treats the 
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FIG. 14: The semiinclusive fourfold differential cross section 
for 3 He(7, pn)p reaction at 55 MeV and 85 MeV photon lab 
energy as function of the np opening angle 8 P + 9 n with 9 V — 
81°. Curves as in Fig. The experimental data arc from 
Ref. El. 
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FIG. 16: 3 He inclusive longitudinal and transverse response 
functions Rl and Rt for the momentum transfer |Q| = 
300 MeV as functions of the energy transfer Op. Curves as in 
Fig. H2 The experimental data are from Ref. |33| (circles) and 
from Ref. [34| (squares). 
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FIG. 15: The semiinclusive threefold differential cross section 
and photon analyzing power E for 3 He(7, n)pp reaction at 15 
MeV photon lab energy as function of the neutron energy E n 
for the neutron scattering angle 8„ = 90°. Curves as in Fig. [(J 



FIG. 17: 3 He inclusive transverse response function Rt near 
threshold as function of the excitation energy E x , Qt being 
the value of three-momentum transfer at threshold. Curves 
as in Fig. El The experimental data are from Ref. |35|. 



screened Coulomb potential between the protons with- 
out partial-wave expansion and therefore has no problems 
with the slow convergence of that expansion. Our solu- 
tion of three-nucleon equations proceeds in partial-wave 
basis and therefore faces the slow partial-wave conver- 
gence of the Coulomb interaction between the charged 
baryons. However, we are able to obtain fully converged 
results by choosing a special form of the screening func- 
tion and by using the perturbation theory of Ref. [l8| 
for treating the screened Coulomb transition matrix in 



high partial waves. This would not be possible, if we had 
used Yukawa screening as in Ref. [6j for two reasons: (a) 
The convergence with respect to screening would require 
much larger radii R; (b) The larger values of R would 
necessitate the solution of the AGS equation with much 
higher angular momentum states. 

(4) Our method for including the Coulomb interaction 
is efficient. Though the number of the isospin triplet 
partial waves to be taken into account is considerably 
higher than in the case without Coulomb, the required 
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computing time increases only by a factor of 3 — 4 for 
each screening radius R, due to the use of perturbation 
theory for high partial waves. 

The obtained results are fully converged with respect 
to the screening and with respect to the quantum number 
cutoffs; they are therefore well checked for their validity. 
The employed technique gets cumbersome in kinematical 
regions with very low relative pp energy, i.e., pp cm. 
energies below 0.1 MeV, due to the need for quite large 
screening radii. 

Physicswise, the Coulomb effect in pd breakup and in 
three-body e.m. disintegration of 3 He is extremely im- 
portant in kinematical regimes close to pp-FSl. There the 
pp repulsion converts the pp-FSI peak obtained in the ab- 
sence of Coulomb into a minimum with zero cross section. 
This significant change of the cross section behavior has 
important consequences in nearby configurations where 
one may observe instead an increase of the cross section 
due to Coulomb. This phenomenon is independent of the 
beam energy and depends solely on specific momentum 
distributions of the three-nucleon final state. Therefore, 
unlike in pd elastic scattering where the Coulomb con- 
tribution decreases with the beam energy until it gets 
confined to the forward direction, in three-body breakup 
large Coulomb effects may always be found in specific 
configurations besides pp-FSl, even at high beam ener- 
gies. 



Another important consequence of this work is that we 
can finally ascertain with greater confidence the quality of 
two- and three-nucleon force models one uses to describe 
pd observables; any disagreement with high quality pd 
data may now be solely attributed to the underlying nu- 
clear interaction. In the framework of the present study 
we reanalyzed the contribution of A-isobar degrees of 
freedom to three-body breakup observables. The largest 
A effects take place in analyzing powers for given config- 
urations. Nevertheless the lack of high quality analyzing 
power data on a broad spectrum of configurations pre- 
vents a full evaluation of the A effects in pd breakup. 
The situation is even worse in three-body photodisin- 
tegration of 3 He, where there are neither kinematically 
complete experiments without polarization nor any ana- 
lyzing power data available. 



Acknowledgments 

The authors thank St. Kistryn and H. Paetz 
gen. Schieck for providing experimental 
data. A.D. is supported by the FCT grant 
SFRH/BPD/14801/2003, A.C.F. in part by the 
FCT grant POCTI/FNU/37280/2001, and P.U.S. in 
part by the DFG grant Sa 247/25. 



E. O. Alt, A. M. Mukhamedzhanov, M. M. Nishonov, 
and A. I. Sattarov, Phys. Rev. C 65, 064613 (2002). 
A. Kievsky, M. Viviani, and S. Rosati, Phys. Rev. C 64, 
024002 (2001). 

C. R. Chen, J. L. Friar, and G. L. Payne, Few-Body Syst. 
31, 13 (2001). 

S. Ishikawa, Few-Body Syst. 32, 229 (2003). 

A. Deltuva, A. C. Fonseca, and P. U. Sauer, Phys. Rev. C 

71, 054005 (2005). 

E. O. Alt and M. Rauh, Few-Body Syst. 17, 121 (1994). 
A. Kievsky, M. Viviani, and S. Rosati, Phys. Rev. C 56, 
2987 (1997). 

V. M. Suslov and B. Vlahovic, Phys. Rev. C 69, 044003 
(2004). 

E. O. Alt, P. Grassberger, and W. Sandhas, Nucl. Phys. 
B2, 167 (1967). 

J. R. Taylor, Nuovo Cimento B23, 313 (1974); M. D. 

Semon and J. R. Taylor, ibid. A26, 48 (1975). 

E. O. Alt, W. Sandhas, and H. Ziegelmann, Phys. Rev. C 

17, 1981 (1978); E. O. Alt and W. Sandhas, ibid. 21, 1733 

(1980). 

R. Machleidt, Phys. Rev. C 63, 024001 (2001). 

A. Deltuva, R. Machleidt, and P. U. Sauer, Phys. Rev. C 

68, 024005 (2003). 

A. Deltuva, A. C. Fonseca, and P. U. Sauer, 
Phys. Rev. Lett. 95, 092301 (2005). 

A. Deltuva, K. Chmielewski, and P. U. Sauer, 
Phys. Rev. C 67, 034001 (2003). 

A. Deltuva, A. C. Fonseca, A. Kievsky, S. Rosati, 



[17 

[is; 

[19 

[20 

[21 
[22 
[23 
[24 
[25' 
[26 
[27 
[28 



[29 
[30 
[31 
[32 

[33 
[34 
[35 



P. U. Sauer, and M. Viviani, Phys. Rev. C 71, 064003 
(2005). 

M. Yamaguchi, H. Kamada, and Y. Koike, 
nucl-th/0310024 

A. Deltuva, K. Chmielewski, and P. U. Sauer, 
Phys. Rev. C 67, 054004 (2003). 

A. Deltuva, L. P. Yuan, J. Adam Jr., A. C. Fonseca, and 

P. U. Sauer, Phys. Rev. C 69, 034004 (2004). 

A. Deltuva, L. P. Yuan, J. Adam Jr., and P. U. Sauer, 

Phys. Rev. C 70, 034004 (2004). 

R. GroBmann et al, Nucl. Phys. 603, 161 (1996). 

G. Rauprich et al, Nucl. Phys. A535, 313 (1991). 

H. Patberg et al, Phys. Rev. C 53, 1497 (1996). 
J. Zejma et al, Phys. Rev. C 55, 42 (1997). 

J. Strate et al, Nucl. Phys. A501, 51 (1989). 

M. Allet et al, Phys. Rev. C 50, 602 (1994). 

M. Allet et al, Few-Body Syst. 20, 27 (1996). 

St. Kistryn et al, Phys. Rev. C 68, 054004 (2003); 

St. Kistryn et al, Phys. Rev. C, to be published, 

micl-th/0508012 

E. Stephan et al, AIP Conf. Proc. 768, 73 (2005). 

R. F. Carlson et al, Lett. Nuovo Cimento 8, 319 (1973). 

W. Tornow et al, AIP Conf. Proc. 768, 138 (2005). 

N. R. Kolb, P. N. Dezendorf, M. K. Brussel, B. B. Ritchie, 

and J. H. Smith, Phys. Rev. C 44, 37 (1991). 

K. Dow et al, Phys. Rev. Lett. 61, 1706 (1988). 

C. Marchand et al, Phys. Lett. B 153, 29 (1985). 

R. S. Hicks et al, Phys. Rev. C 67, 064004 (2003). 



